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\ We study predictions for galaxy cluster observables that can test the statistics of dark 

matter halo shapes expected in a flat ACDM universe. We present a simple analytical 
model for the prediction of cluster-scale X-ray observations, approximating clusters 
as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids 
Q ■ with uniform axial ratios (homeoidal ellipsoids). We test the model against high- 

s' | resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this 

c/3 , simple prescription does a good job of predicting cluster X-ray ellipticities compared 

to the simulations as long as one focuses on cluster regions that are less sensitive to 
recent mergers. Based on this simple model, the distribution of cluster-size halo shapes 
expected in the concordance ACDM cosmology implies an X-ray ellipticity distribution 
' with a mean (ex) = 0.32 ± 0.01, and a scatter a t = 0.14 ± 0.01 for the mass range 

H | (1 — 4) x 10 14 7i M©. We find it important to include the mass dependence of halo 

shape when making comparisons to observational samples that contain many, very 
massive clusters. We analyse the systematics of four observational samples of cluster 
ellipticities and find that our results are statistically compatible with observations. In 
particular, we find remarkably good agreement between two recent ROSAT samples 
and ACDM predictions that do not include gas cooling. We also test how well our 
analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is 
less successful although still useful; the model does not perform as well as a function 
of flux level in this case because of the changing triaxiality of dark matter haloes as a 
function of radial distance. Both this effect and the changing alignment of isodensity 
shells of dark matter haloes leave an imprint on cluster gas that appears to be seen 
in observational data. Thus, dark matter haloes cannot be accurately characterized as 
homeoidal ellipsoids for all comparisons. 
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1 INTRODUCTION 

Clusters of galaxies are the largest bound structures in the 
universe, and the most recently formed ones according to 
the very successful cold dark matter (CDM) cosmology. As 
such, their DM haloes are expected to be less evolved and 
more aspherical than, say, galaxy-size haloes. Most gas in 



clusters DM haloes has not had time to cool, and since it is 
gravitationally subdominant, we can expect it to reflect the 
underlying 3D shape of their dark matter haloes. Indeed, 
large samples of X-ray clusters have been known to show a 
broad distribution of elliptici ties in their surface brightness 
(SB) m aps since the work of iMcMillan. Kowalski. fe UlmeJ 
(1989). A comparison of theoretical predictions to such ob- 
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servations, now that basic parameters of the underlying cos- 
mology are known at the 10%-level or better, may shed light 
on the basic description of the gas in clusters of galaxies. 

The general expectation that in CDM-based theories 
DM haloes are flattened, approximately ellipsoidal, and have 
short-to-long axial ratios as small as s = c/a ~ 0.5 has been 
known for more than 15 years now fearnes fc Efstathioul 
Il987l: iFrenk et aljUiil) . Any asphericity in the DM distri- 
bution has important effects on a variety of observed quan- 
tities. In clusters in particular, asphericity in the dark halo 
potential will map directly to asphericity in the gas density, 
and thus affect the shape of cluster X-ray isophotes and 
Sunyaev-Zel'dovich (SZ) isodecrement contours. Much sub- 
sequent work since these pioneering studies aimed at under- 
standing the origin of such spatial anisotropy, the influence 
of the cosmological model on axial ratios, and improving 
the resolu tion with which the formati o n of DM haloes was 
followed (| pjj^iriskj^_Carlb£rd |l9_9ll; tffiarrgnj ~ 
Ting et alJll99,4 iThomas et al.lll998HSuwa et ".112003 ) . Re- 
cently, higher resolution dissipationless simulations have 
made it possible to fully characterize the scatter and mean 
of axial ratios , as a function of both mass and epoch 
jBullockl 120(4 l Jin g , fc Suto * 120021: iKasun fc Evrardl 12004 
lAllgood et al.ll2005h . I.Ting fc Sutol l|200l hereafter JS) cal- 
culated for the first time axial ratios for isodensity shells, 
using cosmological simulations with 512 3 particles. They 
confirmed that haloes are approximately ellipsoidal in iso- 
density contours, and have provided fits for the dependence 
of axial ratios on mass and epoch. Because the isodensity 
contour method requires a large number of particles, JS re- 
stricted their analysis to haloes more m assive than 6.2 x 
10 12 h ^Mp, in their ACDM simulation. IKasun fc Evrardl 
( 2004) have obtained better statistics for haloes more mas- 
sive than 3 x 10 14 ft^M^ fr om a Hubble volume simula- 
tion. We (Allg ood et al.ll2005l hereafter Paper I) have stud- 
ied haloes spanning the mass range 4x 10 11 — 2 x 1O 14 /i _1 M0 
using several simulations to properly resolve and adequately 
sample the halo population in this entire mass range. Our 
results are consistent with those of JS for haloes of low-mass 
clusters, but yield a steeper mass dependence of axial ratios 
than a simple extrapolation of the scaling relations found by 
JS. This difference is impor tant in the interpret ation of ob- 
servations on galaxy scales llAllgood et aljEo05h and, as we 
show here (Section[3Jl, in the interpretation of X-ray elliptic- 
ities of samples containing ver y massive clusters . Our r esults 
are in agreement with those of lKasun fc Evrardl |2004) if ax- 
ial ratios are calculated in the same manner. However, we 
find here that axial ratios calculated that way are not useful 
for predictions of cluster observables such as X-ray or SZ 
maps (see Appendix^ . A source of uncertainty in the cur- 
rent understanding of halo shapes is the magnitude of the 
effect of gas cooling on clu ster DM haloes iKazantzidis et all 
l2004l:ISpringel et aljEool . 

The variety of current and future observational probes 
of halo ellipticity (see Paper I for a discussion) highlights 
the need to connect these predictions to observations in a 
robust fashion. Here we present an analytic method for pre- 
dicting gravitational potentials and cluster gas density based 
on axial ratios of dark matter haloes. We test the model 
against a sample of 8 high-resolution hydrodynamic ART 
simulations of 7 clusters (mass (1 — 2) x 10 14 /i -1 Mq) and 
1 group (mass 7 x 1O 13 /i _1 M ) in the ACDM cosmology 



llKazantzidis et al"l l2*004'). whose highly variant morphology 
dependence on the line of sight we exploit to statistically 
test the model, and 2 additional clusters from earlier high- 
resolution simulations llKravtsov. Klypin. fc Hoffmanll2002l : 
iNagai fc Kravtsovl koOS). We then apply the model to the 
prediction of ellipticities for cluster X-ray an d SZ maps. 
A pert urbative model has bee n developed by iLee fc Sutol 
( 2003), and further extended bv lWang fc Fan! <2004h to pre- 
dict observed distributions from halo shape distributions, 
but it is not useful for our purposes here because the ellip- 
ticities can be quite large. 

In the following section we discuss the analytic model 
we use to predict cluster X-ray and SZ morphologies, and 
we test it by comparing predictions for morphology of the 
simulation clusters (based on their DM haloes only) with the 
same observable computed directly from the gas density grid 
of the hydrodynamic simulations. In Sectional we focus on 
observations and compare our predictions to several observa- 
tional samples of c luster X-ray ellipticities. We also discuss 
recen t papers (e.g. iFloor et alJl2003l: iFloor. Melott. fc Motll 
l2004h that have compared observed cluster shapes measured 
using X-rays and galaxy distributions to hydrodynamic sim- 
ulations. We finish with a summary of our conclusions. The 
details of the comparison techniques used in Section arc 
outlined in three Appendices:^ Gas Density Inside Triaxial 
Halos; [B] Analytic Potential of Triaxial Generalized NFW 
Halos; and|C| A Comparison of X-ray Ellipticity Measures. 



2 COMPARISON OF MODEL AND 
SIMULATION STATISTICS 

In this section we analyse the prediction for two statis- 
tics of cluster morphology, the mean and the dispersion 
of their ellipticity distribution, expected in a flat ACDM 
universe with Q. m = 0.3, h — 0.7, and o$ = 0.9. We 
first discuss the method to predict cluster ellipticities based 
on their dark matter haloes, and then present the com- 
parison of the predictions to the results from the out- 
put of several high-reso lution hydrodynamical simulations 
(Kazantzidis et al J2004h . In what follows, the virial radius is 
defined as the radius, r v i r , within which the the mean over- 
density drops to A = 18t t 2 + 82(» m (z) - 1) - 39(» m (z) - l) 2 
fervan fc Normanl |l9 98) . Masses are defined as the mass 
within r v i r . 



2.1 Method 

We use a diagonalized moment of inertia tensor iteratively 
calculated within ellipsoids, or ellips oidal shells, to d efine 
axial ratios for dark matter haloes llDubinski fc Carlberd 
Il99ll) . Axial ratios s = c/a and q = b/a (s < q < 1) are 
calculated by diagonalizing the tensor 



^ X i Xj 



R = ^/x 2 + y 2 /q 2 + z 2 /s 2 . 



(1) 



thereby determining s and q for the next iteration. The sum 
is over all particles within a given shell [R, R+AR] , or the el- 
lipsoid interior to R, and the iteration starts with s = q = 1. 
In Paper I we have found that this method predicts axial ra- 
tions that agree with the results of JS, which are based on 
isodensity shells, for cluster-size haloes, provided that the 
axial ratios be calculated within an ellipsoid of semi-major 
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axis R = 0.3r v i r . Here we also find that the same axial ra- 
tios can be used to predict fairly accurately the mean and 
dispersion of the expected X-ray ellipticities, even though as 
often as half the time the predicted ellipticity of an individ- 
ual cluster is off by more than 20%. 

The X-ray surface brightness of an isothermal cluster is 
given by an integral over the gas density squared along the 
line-of-sight (LOS) of a cluster, SB oc J pt as yT oc J p 2 gas . As 
discussed in Appendix [X] under the assumptions of hydro- 
static equilibrium and isothermality, the gas density at any 
point inside a triaxial homeoidal halo can be written in terms 
of the temperature T, the central gas density, and the halo 
potential at the desired point (Eg. lAli . If we assume that 
the halo potential is dominated by the dark matter, then the 
relation is simplified by the fact that the potential of any tri- 
axial generalized NFW halo is analytic ( Appendix^ . Thus, 
using the relation 1A1L we can estimate the X-ray elliptici- 
ties implied by a dark matter halo given the halo axial ratios, 
and an orientation of a LOS corresponding to each of the 
simulation-box axes. 

The only (slight) ambiguity in relating the gas density 
to the halo potential is the the factor T in equation IA1I 
which relates the gas density to the potential exponentially: 
Pgas oc exp(— r$). In Appendix|nj we find that the analytic 
model works relatively well with F ~ c v i r when we compare 
to two high-resolution clusters. This is roughly expected for 
an NFW halo since 



sqpsR] 



. Gfirrip 
kT 



a~ 2 GM vil 



■ /(Cvir 



(2) 



where f(x) = ln(l + x) — x/{l + x) and a is the LOS ve- 
locity dispersion. In the second step we have used sqp s R 2 ~ 
Psphr 2 , where p sp h and r s are spherical-NFW-fit parame- 
ters for the halo, and assumed the expected energy scaling, 
kT ~ nm D a 2 , which is even seen observationally (see e.g. 
iRosati. Borgani. fe Normanll 20021 and references therein), 
albeit with a fair amount of scatter. The final step follows 
from rough scaling relations, and works in detail for the 
clusters we consider in Appendix |C| 1 Therefore, for our 
comparisons we assume 

GM vil fim p 
TVir f(c vil )kT ~ ' W 

and use F = c v i r . For the dark matter halo of the SCDM 
(ACDM) cluster, we find a value of 1.06 (0.98) for the RHS 
of this equation, using the average temperature of the gas 
inside a radius of 400 hT x kpc. 



2.2 Results 

For a given dark matter halo, the method discussed above 
allows us to compute the SB map expected for a given LOS 

1 For example, for the SCDM (ACDM) cluster discussed there, 
G A-f v j r /r v j r /(c v ; r ) = 1063 (630) km/s. The dispersion inside the 
relevant projected radius (400 h~ 1 kpc) for these clusters is sim- 
ilar to these values. For the SCDM cluster, the dispersion is 
a = 1116 km/s (1077 km/s, 1000 km/s) along the x-axis (y- 
axis, z-axis). For the ACDM cluster, we find a = 661 km/s (928 
km/s, 650 km/s). The higher a along the y— axis is due to a merger 
nearly along this axis; however, the same T chosen to fit the ra- 
dial fall-off in surface brightness in the plane perpendicular to 
the x— axis works well for the other two axes. 



through that halo. We discuss in Appendix|n|how an X-ray 
ellipticity can the be obtained from the SB map. There is 
no unique method to calculate ellipticities and, as we dis- 
cuss below and in Appendix|C| it is important to follow the 
procedure chosen by observers to calculate ellipticities in or- 
der to compare to observations. Individual ellipticities can 
differ substantially depending on what part of a map the 
procedure selects and, as we show below, even the means (of 
samples of ellipticities calculated with different procedures) 
will differ. 

Figure^shows SB maps for 8 high-resolution adiabatic 
hydrody namic simulations of clus ters in the ACDM cosmol- 
ogy (see iKazantzidis et alJl2004ll . Each row in the figure 
shows the 3 SB maps corresponding to a LOS parallel to 
each of the coordinate axes of the simulation box containing 
the cluster. We calculate the SB for a given "pixel" in each 
box by summing p 2 gaa \/T over all cells along the LOS-axis. 
Each cell is 7.8 ft -1 kpc on the side, and each map covers 
2 h~ Mpc on the side. The X-ray ellipticity, ex, shown in 
the upper right corner of each map, is calculated using the 
pixels (shown by the shaded areas) containing 20% of the 
total flux above a threshold flux that is 1% of the peak flux 
of the map. This is one of the procedures we consider in this 
work to calculate ellipticities from a SB map. It is a method 
that in the absence of noise yields ellipticities that reflect the 
potential of the DM h alo, as we show in the next paragraph. 
iMcMillan et alJ (jl989fl used this method in their study of X- 
ray ellipticities. However, we find below (see Section|^J that 
their data are heavily affected by noise and do not serve as 
a test of the ACDM cosmology. A more detailed discussion 
of methodologies is presented in Appendix Icl 

For each cluster halo we can use the method described 
in Section f2.1l to compute the predicted SB map for a given 
LOS through a given cluster. We can then compute the pre- 
dicted X-ray ellipticity in exactly the same manner we com- 
puted the ellipticity for the hydrodynamic simulation maps 
of Figure^ In Figure|2]we show a comparison of the two el- 
lipticities. Each point in the figure plots the ellipticity shown 
in Figure Q against the ellipticity calculated based only on 
the DM halo parameters. A Kolmogorov-Smirnov (KS) test 
gives high probability (Pks = 99%) that the two sets repre- 
sent the same distribution. We treat the value of ex for each 
LOS as an independent measurement because for a given 
axial ratio s there is quite a degree of variability expected 
for the other axial ratio, and then there is the variation 
introduced by the orientation of the cluster to the LOS. 
There is indeed quite a degree of variation in ex for each 
cluster, but a stronger test will only be possible with many 
more simulations. The means of the sets indeed agree quite 
well: (e h x ydro ) = 0.36, whereas { e ™° del ) = 0.35. However, 
the dispersions differ significantly: a^ v = 0.18, whereas 

a rnodel _ q rpj^ j g ^ ^ Q ^ middle map m tng lagt 

row of Figure whose ellipticity is greatly enhanced by a 
secondary lump that in this case has a significant relative 
weight due to fact that the flux levels select a narrow region 



0.34, 



of the main cluster. Without that map, {e^ ro ) 
and a^ ydro = 0.15 (recalculating the ellipticity without the 



lump, 



hydro 



) — 0.35 and a 



hydro 



ing difference between a 



hydro 



— 0.15). The remain- 



0.15 and <£' 



= 0.12 



seems to be due to the fact that the gas reflects the chang- 
ing triaxiality in the inner region of DM haloes. We tested 
this by recalculating e ™ odei using axial ratios for the DM 
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Figure 1. X-ray SB maps for hydrodynamic simulation clusters. Each row presents the SB for a LOS along each of the axes of the 
simulation box. Each square is 2 h~ 1 Mpc across. The solid lines show contours of constant SB, spaced by factors of 10. The shaded 
area is the region used to calculate the ellipticity shown in the upper right corners. See text for explanation and discussion. 
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0.2 0.4 0.6 0.8 1 



(F^S f u S 0.01 F p „ k ) 

Figure 2. Model X-ray ellipticities compared to the ellipticities 
calculated directly from SB maps for the hydrodynamic simula- 
tion clusters. Each symbol identifies a cluster, and for each clus- 
ter the symbol plots the ellipticity calculated from the simulation 
map (one for each LOS along each of the coordinate axes) against 
the predicted ellipticity using the method described in Scction l2.ll 
Within the dotted lines the ellipticities differ by less than 20%. 
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Figure 3. Comparison of ellipticities calculated directly from 
SB maps of the hydrodynamic simulation clusters using different 
strategies. The abscissa is the same as in Figure|5] The ordinate is 
an ellipticity calculated using all pixels above a given flux thresh- 
old, in this case 1% of the peak flux in a map. Within the dotted 
lines the ellipticities differ by less than 20%. See text for further 
discussion. 



haloes calculated within R = 0.15r v j r , in which case we find 
(eT del ) = 0.35 and oT° del = 0.15. 

Thus, the analytic model can be used to make fairly 
robust predictions of average X-ray ellipticities. An equally 
robust prediction of the expected scatter does not seem pos- 
sible with a simple homeoidal model of DM haloes, but its 
reliability might be checked by calculating two sets of el- 
lipticities based on DM halo axial ratios calculated within 
two different radii. We have used the model, then, to calcu- 
late the mean and dispersion of X-ray ellipticities expected 
in a ACDM universe at the present epoch. We use a sam- 
ple of 46 DM haloes extracted from the 120 Mpc dis- 
sipationless cosmological simulation discussed in Paper I. 
The cosmology is a flat ACDM universe with £l m = 0.3, 
h = 0.7, and as = 0.9 and the simulation followed 512 3 
particles of mass 1.1 x 10 9 fo _1 Af0. All haloes with virial 
mass (1-4) X 1O 14 ^ _1 M were selected. We calculate their 
axial ratios and concentrations in order to predict X-ray el- 
lipticities for a LOS corresponding to each of the coordinate 
axes of the box. The ellipticity is computed as described 
above, using two flux levels. The samples corresponding to 
each LOS agree quite well with each other. For the combined 
sample wc find 

{e x ) = 0.323 ± 0.013 ; <r e = 0.138 ± 0.008 , (4) 

where the errors are calculated by bootstrap resampling. 
These results are consistent with those for the hydrodynamic 
simulation clusters, for which (ex) = 0.338T0.032 and <r e = 
0.148 ±0.030. 

The strategy to extract an X-ray ellipticity from a SB 
map is by no means unique, and in Figure Figure 0] we 
present two other cases of interest here. For example, in Sec- 
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Figure 4. Comparison of ellipticities calculated directly from 
SB maps of the hydrodynamic simulation clusters using different 
strategies. The abscissa is the same as in Figure |5] The ordinate 
is an ellipticity calculated iteratively using all pixels within an 
elliptical aperture of fixed semi-major axis a = 0.3r v j r . Within 
the dotted lines the ellipticities differ by less than 20%. See text 
for further discussion. 



The Shape of Galaxy Cluster Haloes 7 




Figure 5. Sunyaev— ZeFdovich decrement maps for a hydrodynamic simulation cluster. Each panel presents the map for a LOS along each 
of the axes of the simulation box. Each square is 2 h~ 1 Mpc across. The solid lines show contours of constant temperature decrement, 
spaced by factors of 3. The shaded area is the region used to calculate the ellipticity shown in the upper right corners. See text for 
explanation and discussion. 



tionEJwc discuss a sampl e of X-ray ellipticities obtained by 
iKolokotronis et al.l fcOOll) who use a strategy that empha- 
sizes the central region of a cluster (they were interested in 
mergers). In Figure [3] we show a comparison of ellipticities 
(all calculated directly from the hydrodynamic simulation 
maps) using two different strate gies. The ordinate i s an e l- 
lipticity very similar to that of IKolokotronis et alJ i200ll) . 
calculated using all pixels above a flux threshold correspond- 
ing to 1% of the peak flux of the SB map. The abscissa is as 
in Figure [5] It can be seen there that they differ systemati- 
cally from one another: the means differ by 14%. Therefore, 
a direct comparison of a sample of ellipticities calculated in 
this fashion to our predictions, (Eq. HJ, is not possible. 

Another example of inter e st her e is the strategy used 
bv lBuote. Hart, fc Humphrey! i2005t) . They calculate ellip- 
ticities using all pixels inside a smooth boundary (i.e. the 
boundary is not determined by flux level), which is deter- 
mined by applying the method of lCarter fc Metcalfe! (ll980T) 
(used with out iteration in the stud ies of [McMilla n et alJ 
Il989l and IKolokotronis et"al] l200ll as explained in Ap- 
pendix |nj iteratively, starting from a circle, until the el- 
lipticity converges with a given accuracy. The semi-major 
axis is kept fixed. Figure 0] shows a comparison of elliptici- 
ties (all calculated directly from the hydrodyn amic simula- 
tion m aps) , calculate d using the methodolog y of feuote et alJ 
J2005I) (ordinate) and lMcMillan et alJ (119891) (abscissa). The 
ellipticities agree quite well in mean value and dispersion 
(Pks = 89%), de spite the fact that the methodology of 
iBuote et alJ i2005l) uses all pixels within the elliptical win- 
dow. Thus, the choice of a smooth boundary (rather than a 
flux-selected boundary) makes the ellipticity samples differ 
in no systematic way, unlike the case of Figure |H] A com- 
parison of a sample of ellipticities calculated this way to our 
predictions, (Eq. [IJ, is therefore possible. 

Finally, we also explore here the reliability of the ana- 
lytic model to predict the expected ellipticity of millimeter- 
wave maps of the S Z effect (SZE) in clusters (see e.g. 
ICarlstrom et al]l2002l) . which map the effective temperature 
decrement of the microwave background due to the hot elec- 



tron gas iSunvaev fc Zerdovichll970l) . Figure|S|shows decre- 
ment maps for one of the clusters in Figure (CL5), with 
contours spaced by a factor of 3 (the maps for all clusters are 
included as supplementary material). The maps are quali- 
tatively similar to the SB maps, but the effect of changing 
triaxiality of the DM haloes in the region spanned by the 
isodecrement contours shown is more readily noticed (be- 
cause the signal is proportional p gas instead of Pg as ). We 
show in the upper right corner of each map the ellipticity 
obtained in the same manner as Figure but the decre- 
ment threshold and the percentage of signal in the pixels 
are chosen so that the pixels used cover a region of similar 
size to the corresponding region in Figure^ Specifically, the 
decrement threshold chosen is 10% of the peak decrement in 
the map (as opposed to the 1% of peak signal in Figure 0, 
and the signal in all of the pixels used is 30% of the total 
signal above the threshold (as opposed to the 20% in Fig- 
ure 0. The mean ellipticity and the scatter for this set are 
{(■Sze) = 0.307 ± 0.035 and g S ze = 0.171 ± 0.019, respec- 
tively (bootstrap resampling errors). For the set of Figure 
but with 10% threshold and 30% flux, (ex) = 0.359 ± 0.036 
and ax = 0.175T 0.033, respectively (the difference with our 
1% threshold and 20% flux prediction above is due to the 
changing triaxiality of DM halos). By contrast, the analytic 
model would predict nearly identical destributions. Thus, 
although not as successful as for X-ray ellipticities, the ana- 
lytic model would still be useful to predict, e.g. quantitative 
trends for ellipticities as a function of cluster redshift. 

We now turn to a quantitative comparison of the ex- 
pected distribution of X-ray ellipticities in the ACDM cos- 
mology, obtained using the analytic model we have discussed 
and tested here, with observational samples. 



3 COMPARISON TO CLUSTER-SCALE 
OBSERVATIONS 

Here we compare our predictions to ellipticity distributions 
from samples of cluster X-ray observations. We first anal- 
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yse the methodology employe d by iMcMillan et al] (I1989T) 
and iKolokotronis et al] (1200 ll) to calculate ellipticities for 
their samp les of Abell c luster s. We also consi der briefly the 
sampl e of IMohr et al] (Jl995) considered by IWang fc Far] 
( 2004) for their comparison to observations. These sam- 
ples use different methodology to calculate ellipticities, and 
are affected differently by resolution and noise. Applying 
a KS test to pairs of samples (all converted to 2D ax- 
ial ratios) we find t hat P ks = 0.0031 (Pks = 0-21) for 
the IMcMillan et al. ]l989l) an dlMohr et al] dl995l) samples 
]KolokotrOTuTelTaT (120011) and lMohr et al.l (Il995t) samples') . 
This complicates the comparison of theoretic al predictions 
and o b servations, but it is usually ignored fe.g. lMelott et all 
I 2001] IWang fc Fanl \2004\ . Finally, we analyse a very re- 
cent data set from a nearly co mplete, flux- l imited sample of 
RO SAT clusters discu ssed bv lBuote et al] ]2005t) . 

IMohr et al.l (Il995tl considered a sample of 51 (mostly 
Abell) clusters observed by the Einstein IPC, for which 
they obtained a mean 2D axial ratio, rj, of = 0.80 
and a dispersion a v = 0.12. Converting their axial ratios 
to ellipticities, e = 1 — rj 1 , we obtain (e) = 0.358 ± 0.026 
and <r s = 0.182 ± 0.017. The mean and scatter differ by 
about 1.5 and 2.5 standard deviations, respec tively, from 
our pr edictions (Eq.[IJ. However, the method of lMohr et al] 
( 1995) uses all pixels above a S/N level, and therefore gives 
substantially more weight to the central regions of a SB 
map, where mergers can significantly affect the ellipticity. 
Given our discussion of the results presented in Figure [31 
the difference in mean ellipticity (10%) is entirely within 
the expectection given the different strategy. The agreement 
is somewhat surprising, however, given the potential effect 
that cooling within clus t ers could have on th e DM haloes 
jKazantzidis et al-lliool ISpringel et"ai]|2004l) . We discuss 
this furt her below. 

The IKolokotronis et al. (2001) sample consists of 22 
ROSAT clusters, with a range of velocity dispersions of 
400 — 1000 km/s. Converting their ellipticities to e = 1 — rj 2 , 
the mean and dispersion of their sample are (e) = 0.458 ± 
0.051 and <r E = 0.237 ±0.023. The poor agreement with our 
prediction (Eq. is not surprising given that their method 
emphasizes the cluster centre and there are three clusters 
in the observational sample showing strong evidence of an 
ongoing merger: A2804, A 2933, and A3128 are all bimodal 
iKolokotronis et al-lEoOlh . We have tested that this is in- 
deed the problem by computing ellipticities for the sample 
of hydro clus ters discussed in Sect ion|2] following exactly the 
procedure of IKolokotronis et al] (|S)01), which first defines 
a flux threshold equal to the mean flux within a 600 h -1 
kpc radius, and then uses all pixels above the threshold. A 
KS test between the hydro sample of ellipticities calc ulated 
this way, and the sample of lKolokotronis et al] (1200 1]) . gives 
Pks = 0.82 (the hydro sample is slightly rounder on aver- 
age). Thus we conclude that their sample is in agreement 
with the expectations for a ACDM universe. 

We have also made a comparison with the Einstein data 
of IMcMillan et all (^989) consisting of 49 clusters. Here we 
can expect the comparison to be a better test on the cosmo- 
logical sample because they explicitly exclude image centres, 
thus their shape statistic is less sensitive to mergers (see Ap- 
pendix [CJl. However, the mean and dispersion of their sam- 
ple are (e) = 0.240 ± 0.020 and cr e = 0.142 ± 0.015. There 
is poor agreement with our prediction (Eq. for the mean 
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Figure 6. Simulated X-ray map of a hydrodynamic simulation 
cluster (CL7) after smoothing and adding noise. See text for dis- 
cussion. 



this time. It appears unlikely that this discrepancy could be 
entirely due to missing physics (e.g. cooling) in the simula- 
tions we have used to test the analytic model described in 
Appendix lAl We note t hat even after excl u ding the 3 bi- 
modal clusters from th^jK^lokotroms etal] (|200 ll) sample, 
a KS test against the lMcMlfl£oi^t^dTi^89l) sample (once 
ellipticities are converted to the same definition in terms of 
flux-moment eigenvalues; see Appendix^} rejects that they 
are compatible at the 96% CL. 

The coarser angular resolution of the Einstein data 
probably contributes to this disagreement. For example, we 
find that if we smooth the X-ray map of the ACDM clus- 
ter we discuss in Appendix IU1 (see Figure IUT1 with a gaus- 
sian window of 80 h~ kpc (F WHM, roughly corres pond- 
ing to the 1.6' resolution of the IMcM illan et al.l l|l989l i data 
at the median redshift of their sample, z — 0.057), the X- 
ray ellipticities can change significantly: the entry in col- 
umn 6 of Table E3 would be 0.40 (0.23, 0.46) for the x-axis 
(y-axis, z-axis) as compared to 0.43 (0.28, 0.51) without 
smoothing. This is also co nsistent with the changes found by 
iBuote fc Canizarea ( 1996) for 5 clusters with Einstein, later 
analysed with ROSAT data. We find that a 20% change 
in the predicted ellipticities would make them marginally 
compatible with the data. 

A more important contribution to the difference with 
our p redictions seems t o be t he effect of noise. Many clusters 
in the lMcMillan et al] lll989F) sample have very small elliptic- 
ities but do not look round at all. Figure shows the effect 
of smoothing and noise on one of the clusters in Figure 
(CL7, left). The noise level is at 70% of the threshold (1% 
of peak flu x) used for the calculat ion of ex ■ The automated 
software of IMcMillan et all dl989l) allowed up to 50 islands 
in the X-ray map, thus this map would be accepted. The 
ellipticity without smoothing and noise would be ex ~ 0.18 
in this case, but with smoothing and noise (the latter mak- 
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ing the largest difference) ex = 0.08. If we put this level 
of noise in all the hydro clusters, we find that the hydro 
sampl e becomes fully compatible with the Mc Millan et alJ 
( 1989) sample: Pks = 0.87. Thus, it may be that our pre- 
dictions are not incompatible with this data sample. 

Finall y, we discuss the re cent sample of ellipticities ob- 
tained by iBuote et al. (|2fifl5j) for the flux-l i mited sample 
of B.OSAT clusters of lReiprich fc Bohrineerl l)2002l) . As we 
showed in Section l2.2l a direct comparison to our predictions 
(Eq.0J is in this case possible. It is worth emphasizing here 
that an important advantage of these data is that all the el- 
lipticities are calculated within the same apert ure (0.3r v j r ). 
The m ean and dispersion for the sample of IBuote et alJ 
fcOOifl are (ex) = 0.376 ± 0.019 and <r e = 0.122 ± 0.014, 
respectively (e = 1 — n 2 ). The dispersion is less than ex- 
pected, but only by approximately la. However, the mean is 
substantially larger than expected (by approximately 2.8a). 
This appears to be due to the expected mass dependence of 
axial ratios. 

As discussed in Paper I, dark matter haloes are sys- 
tematically more triaxial the larger their mass. A sim- 
ple relation was found that describes this behavior: (s) = 
0.54(M,/M vir ) - 05 (see Paper I). Many of the clusters in 
the sample discussed bv lBuote et al.l (120051) are much more 
massive than the sample of simulation clusters we used to 
make our predictions (Eq.0J. This is to be expected because 
the observational sample was flux limited and, therefore, 
massive clusters are overrepresented (relative to a volume- 
limited sample, which the simulation clusters represent). Us- 
ing the scaling relation above, we find that our prediction for 
the mean ellipticity of a sample of cluste rs with a mass func - 
tion like that of the clusters analysed bv lBuote et alJ (120051) 
would be (e x ) = 0.353 ± 0.013 (instead of Eq. |]J. 2 The 
remaining difference could well be a statistical fluctuation, 
given that the intrinsic dispersion in ellipticities is a t ~ 0.14. 
Thus, we can expect fluctuat ions Q(a e /\^N) = 0.022 for 
a sample of the size of the IBuote et al] ]2005l) sample. 
By contrast, if we assume the milder scaling advocated by 
JS, (s) = 0.54(M»/M vir ) 03 , the predicted mean elliptic- 
ity would be (ex) = 0.312 ± 0.011 instead, which is 3.4cr 
lower than the observations. It is also worth pointing out 
that the data themselves show evidence of mass depen- 
dence, alth ough not at a high l evel of confidence. If we split 
the data of IBuote et al] ]2005l) by mass, for clusters below 
(above) M vil = 1O 15 /i _1 M (ex) = 0.357 ± 0.027 and cr e = 
0.117±0.020 ((ex) = 0.393±0.026 and a t = 0.123±0.019). 
There are about equal number of clusters in each subsample. 
Although the difference in mean value is not highly signifi- 
cant, it is of the magnitude expected (10%) using the scaling 
of Paper I. 

We conclude from the comparison with these four data 
samples that the predictions for cluster X-ray shapes in the 
ACDM cosmology, assuming gas cooling has only a small 
effect on the shape of their dark matter haloes, are in good 
agreement with the data. A more stringent test, however, 



2 We do the calculation by generating a Monte Carlo set, picking 
the observationally estimated mass of a cluster, and using the 
scaling relation above to get a corresponding (s). We then draw 
axial ratios for the cluster using the form of the distribution of s 
and q found in simulations (see Paper I). 



would require a larger sample of clusters and a better quan- 
titative understanding of the effect of cooling. 

We have attempted to estimate quantitatively the effect 
of gas cooling on cluster X-ray ellipticities, which generi- 
cally makes DM haloes less triaxial. We use a hydrodynamic 
simulation of one cluster for which cooling and star forma- 
tion was abruptly terminated at redshift z — 2 in order for 
the cluster to have rea sonable star and gas fractions (see 
iKazantzidis et alJl2004T) . The effect of cooling on DM halo 
axial ratios for this cluster agre es very well with the aver- 
age effect shown in Figure 4 of ISpringel et all ]2004]) . We 
calculated the short /long axial ratio in logarithmic radial 
distance bins i n ord er to directly compare to the figure in 
ISpringel et al] i2004T) . We find that there is good agreement 
in the distance range (0.1 — 0.3)r v i r . Therefore, we have esti- 
mated the expected effect on X-ray ellipticities in two ways. 
We can compute the change in ellipticity by comparing the 
ellipticities with and without cooling for this one cluster. 
Since the change in axial ratios seems to be representative 
of the expected average change, we can estimate that the 
effect would be to make X-ray ellipticities 10-20% smaller. 
We have also estimated the effect by generating a catalog of 
DM "haloes", where a halo is represented as a set of axial 
ratios and a concentration. We generate axial ratios using 
the form of the distribution of s and q found in simulations 
(we use s = 0.54 and a 3 = 0.1 for the gaussian distribution 
of s; Paperl). We g enerate concentrations using the lognor- 
mal distribution of IWechsler et all ]2002l) . with mean of 7 
and a log-dispersion of 0.14. Finally, we orient randomly 
the principal axes in a box. A mean short/long axial ratio 
s — 0.54 instead of s = 0.45 adequately rep resents the ef- 
fect seen on average by ISpringel et all |2004), and the effect 
on the cluster discussed here. We find in this case that we 
can expect X-ray ellipticities to be ~ 25% smaller. It is thus 
rather surprising that we find the level of agreement with 
the data that we have found here without taking the effect 
of cooling into account. 

We have considered a lower-as cosmology, in which DM 
halos are predicted to be more triaxial (see Paper I), as a 
possible explanation of this surprising result. If DM halos 
were more triaxial, the predicted X-ray ellipticities would 
increase and then cooling could bring results into agreement 
with the data. We have found in Paper I that a simple scaling 
relation accounts for the dependence of axial ratios on as 
(see Paper I, Eq.(7)). The predicted s can then be used as 
above to generate a catalog of axial ratios. We find that even 
for a value of as as low as as = 0.75, the expected m ean X- 
ray ellipticity of a sample like the lBuote et al] ]2005l) sample 
changes only to (ex) = 0.378 ± 0.013 (from (ex) = 0.353 ± 
0.013 for as = 0.9). 

There are potential biases that can affect comparisons 
of the model with observations. For example, in relaxed 
cooling flow clusters the temp erature decreases toward the 
center in the cluster core Ce.g.. lDe Grandi fc Molendll2002l : 
Vikhlinin et al .l2005aD . Line emission of low-temperature X- 
ray gas can significantly alter the p^asVT weighting assumed 
in our analysis and, therefore, the shape of the X-ray bright- 
ness isophotes. To take this effect into account, however, we 
need to know the temperature and metallicity distribution 
in clusters. We plan to address this issue in future work, 
which will make use of the cluster simu lations with galaxy 
formation. However, IBuote et al] (120051) have calculated el- 
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lipticites in annuli as well, i.e. excluding the cluster centers 
altogether. The mean ellipticity is only slightly (4%) higher. 

Another possible source of bias can arise in comparisons 
with the shape estimates based on the isophotes defined at a 
constant fraction of the peak flux of the cluster. The profiles 
of real clusters are ofte n quite "cuspy" in their centers (e.g., 
IVikhlinin etafl l2005bl) and are considerably steeper than 
the radial gas density profiles of clusters in our adiabatic 
simulations. This difference in the radial gas distribution 
will result in different radii of the isophotes defined with 
respect to the peak flux. This may mean that the shapes 
would be measured at systematically different radii in sim- 
ulations and obser vations (smal l er rad ius in observations). 
Note, however, that Buote et all i2005f) calculate ellipticities 
within the same apperture, as we have stressed above. 

Two qualitative trends in X-ray maps appear to reflect 
the more complex nat ure of dark matter haloes s een in high 
resolution simulations. iBuote fc Canizaresl lll996t and refer- 
ences therein) have pioneered detailed studies of X-ray maps 
to constrain cluster halo profiles. They studied 5 Abell clus- 
ters using oblate and prolate spheroids in order to bracket 
the possibilities, and concluded that the ellipticity (there 
is only one axial ratio if one assumes oblate and prolate 
spheroids) of the haloes was constrained to be in the range 
0.40 — 0.55. The systematic trend of interest here is that 
4 of the 5 clusters show a decreasing ellipticity of the X- 
ray isophotes at larger radii. A similar trend can be seen 
in the gas data for the two high-resolution simulation clus- 
ters discussed in AppendixlUl (see Figure lUTl and Figure lU2t 
and is due to the decreasing triaxiality of the dark matter 
halo at larger radius. The effect is not very pronounced, so 
the simple isothermal/homeoidal halo model could still be 
used for the ellipticity comparison above. The same is not 
true, however for Sunyaev-Zel'dovich decrement maps (see 
Figure |C3j. As discussed in Appendix IU1 the different sen- 
sitivity of SZ maps to density and temperature (oc J p gaa T) 
make these observations more sensitive to our simplistic as- 
sumptions, and the simple mapping from halo shape param- 
eters will break down more visibly. More detailed modeling 
will likely be required to interpret SZ shape measurements 
accurately. 

A second c o mplic ation of interest here was noted by 
iMcMillan et alJ il989^ . who pointed out that a fraction 
(~ 15%) of their clusters exhibited isophotal twist with 
a "continuous rotation of the intermediate isophotes". We 
have found that one of the high-resolution simulation clus- 
ters discussed in Appendix |C| (see Figure ICJ2ft shows this 
kind of twist due to coherent twist of the dark matter den- 
sity shells. Of course, it will be interesting to quantify the 
frequency of this effect, as well as its origins. The degree 
of misalignment in the case of this cluster (~ 45° in the 
radial range ~ (0.3 — l)r V ir) would be rare judging by the 
results of JS for 12 haloes. However, a direct comparison is 
not possible because the angles involved are not the same. 

Therefore, while some observational quantities are 
somewhat insensitive to the complex non-homeoidal nature 
of halo structure, many observed properties are quite sen- 
sitive to changing ellipticities and twists. Specifically, the 
higher-order trends in halo shapes may leave imprints in 
cluster gas that could be studied in detail by analyses of 
X-ray and Sunyaev-Zel'dovich maps. 

Clusters X-ray ellipticities can be expected to evolve 



with redshift due to increased halo triaxiality (see Paper 
I, and references therein). Recent papers have called atten- 
tion to a possible significant evolution of the ellipticity with 
redshift even over the nearby redshift range z = — 0.1 
jMelott. Chambers, fc Milled! 200ll|Plionisll2002l) . and have 
claimed that cluster X-ray and optical profiles are a little 
less flattened than pred icted by dissipationless and hydrody- 
namic simulations fsee lFloor et al.ll2003l : iFloor et al1l2004l: 
and references therein). However, it is important to compare 
observational data to simulated clusters of similar mass (the 
Floor et al. clusters were more massive than most of the ob- 
served clusters) and, as we have explained (see AppendixICt. 
to mimic the way the data was treated. It is hard to draw 
clear concl usions when the rat h er dif ferent IMcMillan et alJ 
i 19891) and iKolokotronis et alJ <l200ll) X-ray data sets an d 
analyses are combine d, as was done bv lMelott et al 1 1 20011) . 
Ijeltema et al] <l2005l) have detected evolution in cluster mor- 
phology in a more homogeneous sample of clusters with 
Chandra data. However, the morphology is not quantified 
as an ellipticity, therefore we cannot assess how well this 
observation constrains theory in our current paper. 



4 CONCLUSIONS 

We have presented a simple analytic model to predict clus- 
ter halo gas profiles based on dark halo shapes, under the 
assumption that clusters are isothermal and in hydrostatic 
equilibrium within haloes that are homeoidal ellipsoids (i.e. 
with constant axial ratios). We found that certain observa- 
tional properties, such as ellipticities of X-ray surface bright- 
ness maps, can be adequately described by this modeling. 
When applied to our sample of cluster-mass haloes we find 
that the predicted distribution is in good agreement with 
observational samples of el lipticities for galaxy clusters. The 
agreement with the recent lr3uote et aTM2005T) analysis of a 
complete ROSAT sample is especially impressive, because 
we found it important to take into account the predicted 
mass dependence of halo shape in comparing to this data 
sample rich in very massive clusters. The usefulness of our 
model is sensitive to how the observed ellipticity is defined. 
Specifically, care must be taken to avoid definitions that 
make the calculated ellipticities sensitive to mergers. 

The shape of dark matter haloes undoubtedly cannot 
be fully characterized by simple models with constant axis 
ratios. While we have used inertia-tensor derived axial ra- 
tios to characterize halo shapes in a simple way, and ex- 
plored how simple assumptions about halo shapes can be 
used to compare to observational tracers of halo structure, 
we find that more detailed predictions will be required for 
many observational comparisons (Appendix |CJl For exam- 
ple, the isothermal/homeoidal assumption becomes less use- 
ful for comparison to measurements like Sunyaev-Zel'dovich 
decrement maps. In addition, radially decreasing ellipticities 
can arise from the changing shape of isodensity contours 
with radius, and twists in X-ray isophotes can arise from 
misalignment of contours at large and small radius. Predic- 
tions aimed at this kind of data will require a more detailed 
analysis of ACDM halo shapes, including a detailed char- 
acterization of halo ellipticities as a function of radius, and 
the frequency of isophotal twists. Work in this direction is 
under way. 
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APPENDIX A: GAS DENSITY INSIDE 
TRIAXIAL HALOS 

Here we present a simple model of the gas density expected 
inside a cluster halo, and use it to calculate X-ray proper- 
ties such as surface brightness. The model can also be used 
for other gas-density-dependent observations, such as SZE 
maps from millimeter- wave observations of clusters (see e.g. 
ICarlstrom et al]l2002fr . We define the dark matter halo den- 
sity model, and calculate its potential, in Appendix [5] The 
gas density model is based on three common approximations 
about the gas: 

1) the gas is in hydrostatic equilibrium, 

2) the gas is isothermal, 

3) the gas makes a negligible contribution to the total mass. 
These assumptions are quite restrictive, although it is triv- 
ial to modify equation 1AH below for a polytropic gas. In 
Appendix[U]we shall relax all of the assumptions and work 
directly with the gas in two high-resolution simulations of 
galaxy clusters. We work out expected X-ray properties for 
the clusters in the simulations, and compare them to the 
predictions based on the model described here. We find that 
the model works fairly well, despite its simplifying assump- 
tions. We further test the model statistically against a small 
sample of high-resolution simulation clusters in Section |5] 

With the assumptions listed above, the gas density ex- 
pected at a point (x, y, z) inside a triaxial halo can be written 
in terms of the gas density at the origin, the dark matter po- 
tential $(x,y, z), and the gas temperature T. For the halo 
density model discussed in Appendix [B] we find it conve- 
nient to work with the potential in units of the overall factor 
AitGsqpsR 2 . Therefore, we write 

p gas (x,y,z) = p gas (0)exp|-r (<b(x,y,z) - $(0)j} . 

(Al) 

Here &(x,y,z) = &(x,y, z)/iirGsqp s R 2 , so the constant T 
is given by 

T = ^Gs q p s R 2 s & , (A2) 

where p is the mean molecular weight. For clusters 
with galaxy velocity d ispersion a, kT ~ pm p a 2 (see 
e.g. iRosati et al] Qi002). Therefore, since 4irsqp a R 2 = 
0(a 2 /G) c v i r (see Section^, we can expect F ~ c vlv . 

We can use this simple model to calculate the ex- 
pected X-ray SB of hot gas in a dark halo with a given 
potential. Since we have assumed the gas is isothermal, SB 



oc f p' gas , where the integral is calculated along the LOS. 
In Appendix 151 we calculate the potential &(x,y,z) in the 
principal-axis coordinate system of the dark matter halo. 
Therefore, in order to calculate SB we need the orientation 
of the LOS in this co ordinate system. We use the conven- 
tions of iBinnevI il985T> . in which the LOS-axis is defined by 
azimuthal and polar angles <j> and 8, respectively. 

We thus have the following model to predict the X-ray 
SB map expected for a given projection of a dark matter 
halo in a simulation box. We first calculate the axial ra- 
tios s < q < 1 by the iterative procedure described in Sec- 
tion |21 in Appendix |C| we find that axial ratios calculated 
using a solid ellipsoid of semi-major axis 0.5r v i r works well 
to predict flux-weighted ellipticities. We also obtain from 
the procedure the orientation (<f> and 6) of a given LOS, and 
the orientation (position angle, PA) of the projection of the 
shortest axis of the halo on the plane perpendicular to the 
LOS. For a given point along the LOS, we find p gas by first 
rotating its coordinates in the plane by the PA. We then 
apply the inve rse of the rotation parametrized by (f> and 9 
( Binnevl ll985l) . This gives us the coordinates of the point 
along the LOS in the principal-axis system, from which we 
obtain p gas using equation iAll . Therefore, we can calculate 
J p 2 as numerically at any given point on the plane. We shall 
refer to this model for the SB as the "analytic model" (even 
though it involves numerical integration), in order to distin- 
guish its predictions from those we work out directly from 
the gas density in two high-resolution simulations of galaxy 
clusters that we analyse in Appendix [U] where we compare 
predictions for X-ray ellipticities. 



APPENDIX B: ANALYTIC POTENTIAL OF 
TRIAXIAL GENERALIZED NFW HALOS 

Here we consider the potential of triaxial dark matter haloes 
with a density profile that is a simple generaliz ation of a 
specia l case of the spherical profile introduced bv lHernauisj 
(1990). We assume that isodensity shells are homeoidal el- 
lipsoids, i.e. with constant axial ratios s and q (s < q < 1), 
and that the density profile is given by 

/ \ _ Ps 

P[x,y,Z) ( R / Rs ) a ( 1 + R / Rs ) v -a 

R= ^x 2 + y 2 /q 2 + z 2 /s 2 (Bl) 

where p 3 and R s are a scale density and radius, respec- 
tively. Assuming constant axial ratios allows us to reduce 
the calculation of the potential to a one-dimensional inte- 
gral, which in some cases can be solved analytically, after 
a simple transform ation of the general re sult for ellipsoidal 
mass distributions llChand rasckhar 1969). 

We first consider 77 = 3. This was found to be a 
good approximatio n (assuming constant axial ratios) by 
Ijing fc Sutol i2002T) for their 12 high-resolution haloes. We 
also find this to be a good approximation for the haloes of 
two high-resolution hydrodynamical simulations we discuss 
in Appendix|n| However, since spherical fi ts to large samples 
of haloes find deviations f rom this value iAvila-Reese et al] 
ll999UThomas et aljllool . we also generalize the results to 
other values below. 
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The potential of a thin homeoid of mass M (axes a > 
b > c), at a point (x,y,z ) outside the shell, can be written 
as l)Chandrasekharl (l969 ) 



<P M (x,y,z) = — 



da 



2 Jx ^(a? + «) ( b2 + u ) ( c2 + u ) 

(B2) 

The parameter A in equation l|B20 is the parameter of the 
confocal ellipse passing through (x, y, z); it is the largest root 
of 



2 2 2 

x y z 

+ ,n ; + — r = 1. 



a 2 + X ' b 2 + X ' c 2 + X *• (B3) 

Since the integral 1B21 can be solved analytically, we find 
that 



$> M (x,y,z) = - 



GM EllipticF f /a2 ~ c2 /a2 - &2 



Va 2 - c 2 ^ V a 2 + A ' V a 2 — c 2 

(B4) 

The potential insid e the homeoid is a constant 
iChandrase khar 1969/) , therefore it is given by <$>m(x, y, z) 
with A = 0. 

We can construct the potential inside a triaxial NFW- 
type halo now, assuming homeoidal symmetry (i.e. constant 
axial ratios). First, for the potential at (x,y,z) due to all 
mass shells inside (i.e. inside the shell passing through the 
point) we find 



$in = -Af- 01 I dm - g-_ 



(1 + mC) 3 - 



EllipticF 



1 - s 2 



l + A(m)/m 2 .R 2 ' V l-s 



1 - q 2 



(B5) 



Here ( = R/R s , A = 4nGsqp 3 Rs/Vl - s 2 , and A(m) is the 
largest root of 



+ ■ 



y 



+ 



= 1 . (B6) 



m 2 R 2 + X m 2 R 2 q 2 + X m 2 R 2 s 2 + A 
Second, for the potential due to all shells outside we find 



<&out 



2- a 



EllipticF ^y/l - s 2 , 



i-r 



i + C 



(B7) 



The total potential is then $(a;, y, z) = "l> m + $ ou t • Also, 
we find that the only change needed for r\ 7^ 3 is to replace 
(l + m() 3 ~ a by (l + m() r '~ a in equation 1B5L and to replace 
(1 - (C/(l + 0) 2_Q ) in equation {Eft by 

(2 - a) T(2 - a) T(r] - 2) 



r(r, - a) 

C 2 -" hypergeom([l, g - a], [3 - a], C/(l + 0) rR8 v 

(i+o— ■ (B8) 

For spherical symmetry (q = s — 1) we can check the 



standard result for the potential of an NFW halo. We can 
use EllipticF(x, 1) = x + Q(x 3 ) to get 



: EllipticF 



l-s 2 



1 - q 2 



1 + X{m)/m 2 R 2 ' V l-s 2 
1 



/I + A(m)/m 2 i? 2 



(B9) 



In this case (q — s = 1) the right hand side is just m, and 
for the NFW profile (a = 1, 77 = 3) we have 



1 mdm _ ln(l + C) + ln(l + C) C - C 



(1 + ™C) 2 



C 2 (l + C) 



(BIO) 



Therefore, at radial distance r — x 2 + y 2 + z 2 
2 (\n(l + r/R s ) 1 



$m(r) = -AnG p a Rs 



r/Rs 1 + r/Rs 



®out{r) = -AnGps R s 



1 r/Rs 



1 + r/R s 

and the total potential takes the standard form, 



(Bll) 



H r) = -4,G Ps R/^l+im. (B12) 
r/Rs 

For a/1 and rj ^ 3, the integral in equation 1B5II can be 
obtained analytically. The spherical potential of a general- 
ized NFW halo is then 



$(r) = -4nG P sR, 



where 



a f r(2 - a) T( V - 2) 

1 r( v -a) 

(r/R s f- a A(r/Rs) 
(2-a)(l + r/R a )i- a 



(B13) 



2 q / 2; 

A{x) = hypergeom I [1, r] — a], [4 — al 



3 — ol 



J ' 1 + x 



hypergeom [1, 77 — a], [3 — a] 



1 + x 



(B14) 



APPENDIX C: A COMPARISON OF X-RAY 
ELLIPTICITY MEASURES 

Here we evaluate the reliability of our method for predicting 
individual cluster X-ray ellipticities based on the calculated 
axial ratios of dark matter haloes. Using a h ydrodynamical 
simula tion (with cooling) of a single cluster iBuote fc Tsail 
( 1995) found that this method, assuming isothermal gas, 
allows an accurate estimation of the ellipticity of the dark 
matter even if the gas has a strong temperature gradient, 
so long as any substructure in the cluster is excluded in the 
analysis. We use high-resolution adiabatic hydrodynamical 
simulations of two clusters in order to calculate X-ray SB 
maps directly from the gas data of the simulations. We then 
compare these maps in detail to predictions based on the 
properties of their dark matter haloes, using the theoretical 
model described in Appendix^ We find that the theoretical 
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Table CI. Ellipticity results for ACDM cluster; see Figure^] Here e x = 1 - A^/A 2 . 



LOS 


gas 

e x 
> 0.01 


model 

> 0.01 


gas 
t X 

> 0.1 


c model 
t X 

> 0.1 


gas 
t X 

0.01 - 0.08 


^model 
t X 

0.01-0.08 


gas 
L X 

0.1 - 0.2 


c model 
t X 

0.1 - 0.2 


x-axis 


0.61 


0.48 


0.79 


0.49 


0.43 


0.46 


0.72 


0.53 


y-axis 


0.24 


0.23 


0.40 


0.24 


0.28 


0.23 


0.32 


0.26 


z— axis 


0.65 


0.51 


0.77 


0.53 


0.51 


0.50 


0.69 


0.48 



model can work relatively well (predicting ellipticities within 
10% of the gas-data values) depending on exactly how the 
observational ellipticity is defined. A statistical (rather than 
case-by-case) test of the model is presented in Section |5] 

We fir st discuss a ACDM cluste r that has been studied 
in detail bv lNaeai fc Kravtsovl (120031) . In Figure lCTl we show 
a "surface brightness" map calculated from the gas data of 
the simulation. The solid lines really show contours of con- 
stant value of J Pj as , where the integration is along a LOS 
parallel to the x-axis of the simulation box. Of course, SB 
K f PgasVT, but we have dropped the temperature depen- 
dence for simplicity, given that it makes only a small differ- 
ence in calculated ellipticities (;$ 5%; see below). We calcu- 
late SB for a given "pixel" by summing over all cells along 
the LOS. Coordinates are shown in pixels, with 7.8 h~ kpc 
per pixel. The dotted-line contours are the SB contours pre- 
dicted by the model described in Appendix[X]with the factor 
r chosen to match the radial SB profile, F = 10.5. For this 
halo Cvir = 11.5, therefore T ~ c v i r , as expected (see Sec- 
tionEJ. TheX -ray ellipticities discussed here are only mildly 
dependent on F (e.g. e% odel = 0.46 in Table IUT1 changes to 
e modei _ o.40(0.48) for V = 8(13)). They are mostly sensitive 
to the axial ratios s and q, and the relative orientation of the 
LOS, described by polar angles (f> and 6 in the principal-axis 
coordinate system. The axial ratios and polar angles were 
calculated inside an ellipsoid of semi-major axis 600 kpc 
using the iterative method described in the text, and are 
shown at the top of the figure. The dashed-line contour il- 
lustrates a predicted isophotal con tour based on axial ratio s 
calculated with the prescription of iKasun fc Evrardl i2004f) . 
The result is very similar for SZ isodecrement contours. 

The highly irregular, innermost solid contour is due to 
a minor merger nearly in the "plane of the sky" (about 25° 
off the y-axis of the box). The merger is ideal to test quan- 
titatively observational strategies to calculate an ellipticity 
that best represents the global triaxiality of the dark matter 
halo. It is also an ideal test of the reliability of our method 
because it allows us to gauge what bias can be introduced 
in the calculation of ellipticities by the presence of a minor 
merger, which can be expected to be common for cluster-size 
systems. We discuss both of these issues below. 

There are various strategies to calculate X-ray 
ellipticities. Here we c onsider the method used by 
Kolokotronis et al.l |200J> (22 clusters; ROSAT data) and 
McMillan et alJ l)l989T) (49 clusters; Einstein data) as exam- 
ples used in anal yses of samples of cluster s. Both studies 
use the method of ICarter fc Metcalfe! il980f) adapted to an 
X-ray image. The ellipticity is calculated from the positive 
roots A+ and A_ (A+ > A_) of the characteristic equation 

(/x 20 -A 2 )(/xo2-A 2 ) = //? x . (CI) 
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Figure CI. Surface brightness plot for the ACDM cluster in the 
yz-plane of the simulation box. The axes are in pixels (7.8 h^ 1 
kpc per pixel) and the solid lines show contours of constant J pg as , 
spaced by factors of 10. The SB peak (centroid) is indicated by the 
cross (open circle), and the innermost solid contour corresponds 
to a level of 15% of peak value. The dotted-line contours show the 
predictions of the analytic model described in Appendix |X] The 
factor r is estimated by fitting, for one projection, the radial SB 
profile of the simulation. It is then used for all other projections. 
Finally, the dashed-line contour illustrates a predicted contour if 
a more global measure of triaxiality were used. See the text for 
further discussion. 

The moments p mn are defined in terms of the flux at a 
given pixel (a*,j/j) by 

ij ij 

where (x,y) is the image centroid (x = x ifa / J2ij fa> 
y = ■ Djfij I ■ fij)- The ellipticity is then calculated as 

e x = l~Al/A 2 + (C3) 
bv lMcMillan et alJ Jl989h . and as 

ex = l-A-/A+ (C4) 
by iKolokotronis et alJ i200ll) . We shall use equation 1C3I 
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here, except when comp aring directly with the data of 
IKolokotronis et al.l (120011) . 

The X-ray ellipticity ex is rather sensitive t o what pix- 
els are used to calculate it. IKolokotronis et al.l i200lf) use 
all pixels above a flux threshold (which is the average flux 
within a region of given radius). For example, in Figure IUT1 
this threshold is ~ 0.01 of the peak flux within 600 h -1 kpc 
(which is the largest radius they use to define the threshold). 
In TableEI] (columns 2-5) we show results for the cluster of 
Figure IUT1 for two flux thresholds (0.01 and 0.1 of the peak 
flux) and for a LOS along each of the axes of the simulation 
box. This choice of pixels emphasizes the brightness peaks 
and, therefore, is more sensitive to mergers. Thus, the an- 
alytic model prediction for the ellipticity, e ™ odei ; deviates 
significantly from the value calculated directly from the gas, 
t 9 x S , except w hen the merger i s near ly along the LOS. On 
the other hand, McMillan et al. ( 1989) explicitly exclude the 
centre of an image in order to characterize the global dynam- 
ics of a cluster. They use all the pixels containing 20% of the 
flux above a faint threshold. The latter varies substantially 
across the sample, but for 80% of clusters it is ~ 0.01 — 0.2 
of the peak flux. In Table lUTl (columns 6-9) we show results 
for two flux ranges (0.01-0.08 and 0.1-0.2) and for a LOS 
along each of the axes. In this case the analytic model per- 
forms much better, provided the fainter threshold is chosen 
low enough. If we calculate the deviation of the model ellip- 
ticity from the gas ellipticity (using a flux threshold of 0.01) 
for 100 rand om LOS we find 

(1) For the iMcMillan et alJ Jl989l) ellipticity the analytic 
model works fairly well; 2/3 of the time the model predicts 
the ellipticity within 10% of the gas value. Also, it would 
not bias a statistical sample because it predicts larger and 
smaller val ues with equal frequency. 

(2) For the lKolokotronis et alJ i200lt) ellipticity the analytic 
model predicts a systematically smaller value. This is ex- 
pected in this case because the model misses the merger, 
therefore it predicts rounder SB contours from all viewing 
angles. In this case we find that 2/3 of the time the value is 
20 - 30% smaller. 

We find similar trends for the mean and the dispersion of 
ellipticities calculated for the sample of clusters discussed in 
Section [5] although individual values can deviate more than 
indicated here. 

The flux level at the outermost contour in Figure IU1I 
is ~ 0.002 of the peak flux. At this flux level the con- 
tour is clearly rounder than the model prediction (due to 
the fact that the dark matter halo gets rounder farther 
out, whereas the analytic model assumes constant axial ra- 
tios). However, pixels up to much higher flux levels (~ 0.06 
of peak flux) enter the calculation in order to accumulate 
20% of the flux above t his fa inter threshold in the ap- 
proach o f lMcMillan et al.l il989T> . For example, for the x-axis 
e 9? s = 0.41 and e^ odel = 0.45 in the flux range 0.002 - 0.06 
of peak flux. Therefore, the model works well down to lower 
thresholds. 

The analytic model assumes that the gas is isothermal 
in order to predict the SB. We can check how much this 
is likely to affect a comparison with actual data by calcu- 
lating the ellipticity from the simulation data including the 
temperature dependence. We find that for the flux levels 
considered here, the effect is rather small. For example, the 
entry in column 2 of Table would be 0.59 (0.23, 0.63) for 



the x-axis (y-axis, z-axis) as compared to 0.61 (0.24, 0.65) 
assuming isothermality. The temperature in this cluster falls 
by a factor ~ 1.9 in the radial rang e (0.1 — 0.5)r v j r , which 
is co nsistent with observations (see iDe Grandi fe Molendil 
I2002L and references therein). Therefore, the temperature 
variation of the simulation gas is representative of that of 
real clusters. 

We have also tested whether using the dark matter po- 
tential of this cluster directly would significantly improve 
the prediction for ex- The assumptions are still the same, 
but the potential is calculated directly from the dark mat- 
ter distribution in order to predict the gas density. We find 
that the results improve as follows. For example, the entry 
in column 3 of Table IUT1 would be 0.53 (0.25, 0.55) for the 
x-axis (y-axis, z-axis) instead of 0.48 (0.23, 0.51). 

Finally, in order to study whether the analytic model 
indeed performs better in the absence of a merger, we have 
analysed in the same manner a high-resolution simulation 
cluster that does not have an ongoing merger . It is a SCDM 
cluster discussed in detail bv lKravtsov et all <|2002T) . In Fig- 
ure lC2l we show the SB map calculated as in Figure IUT1 and 
in Table IU2l we show the results for the ellipticity. For this 
cluster r = 9.3 and c v i r = 10.4, therefore F ~ c v j r as before. 
In this case the model works reasonably well for either one 
of the definitions of ellipticity, provided that the flux thresh- 
old is sufficiently high. For faint thresholds the model fails 
to reproduce the trend of rounder and twisted SB contours 
in the simulation (for LOS = y-axis, the reverse trend in 
Table IC2I is due to the chance projection of a distant hot 
spot that appears only at a level ~ 0.01 of peak flux). It is 
in fact the twisted SB contours that cause most of the differ- 
ence between model and simulation gas. This is due to the 
fact that isodensity shells are fairly misaligned in this case. 
The projected, 100 kpc-thick isodensity shell of 400 h~ x 
kpc (800 - 900 hT x kpc, 1000 - 1050 h' 1 kpc) semi-major 
axis makes a 15° (35°, 60°) angle with the vertical direction 
in Figu re | U2 | Such l arge m isalignments were found to be 
rare bv lJing fc Sutol (I2002T) . therefore we assume here that 
the model also works down to the faint threshold level of 
~ 0.01 of peak flux in the absence of a merger, and for both 
ellipticities. 

In Section|3]we analyse the expected distribution of X- 
ray ellipticities for cluster-mass haloes in the cosmological 
box discussed in Section |5] We calculate ellipticities using 
th e analytic model , and c omp are the distribution to the d ata 
of IMcMillan et al.l Jl988f) and IKolokotronis et all <200llh 

We have also considered the reliability of the analytic 
model to predict the shape of SZE maps of clusters. In Fig- 
ure lC3l we show "temperature decrement" maps for the two 
clusters we have discussed. The solid lines show contours 
of constant value of J p gas T (spaced by factors of 3). The 
integration is along a LOS parallel to the x-axis of the corre- 
sponding simulation box, as was the case in Figures IC II and 
IC2I Since the dependence on gas temperature is linear in 
this case, we can expect a more significant effect of temper- 
ature on the shape of isodecrement contours. The top panels 
of Figure IU31 compare the shape of the contours of constant 
value of J p gas only (dashed lines) to decrement contours. 3 

3 The dashed-lino contours are not shown spaced by a fixed fac- 
tor. The levels are just chosen to give contours of similar size to 
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Table C2. Ellipticity results for SCDM cluster; see FigureEU Here e x = 1 - A?_/Ai 



LOS 


gas 

e x 
> 0.01 


model 

> 0.01 


gas 
t X 

> 0.1 


c model 
t X 

> 0.1 


gas 
t X 

0.01 - 0.09 


,modeZ 
t X 

0.01-0.09 


gas 
L X 

0.1 - 0.2 


c model 
t X 

0.1 - 0.2 


x-axis 


0.30 


0.35 


0.38 


0.38 


0.25 


0.34 


0.38 


0.34 


y-axis 


0.48 


0.37 


0.45 


0.42 


0.51 


0.35 


0.42 


0.39 


z— axis 


0.48 


0.55 


0.52 


0.56 


0.44 


0.54 


0.48 


0.55 




Figure C2. Surface brightness plot for the SCDM cluster in the 
yz-plane of the simulation box. The axes are in pixels (15.6 h 
kpc per pixel) and the meaning of symbols is as in Figure IC1I 
In this case the innermost solid contour corresponds to a level of 
25% of peak value, and each solid contour is drawn at one tenth 
of the solid— contour levels of Figure |CT1 

Both sets of solid contours are calculated directly from the 
gas and temperature data of the corresponding simulation. 
It can be seen that in the presence of a merger (the ACDM 
cluster case) the temperature dependence indeed makes the 
isodecrement contours noticeably different from contours of 
J Pgas- However, in the absence of a merger (the SCDM clus- 
ter case) they agree fairly well in shape. For example, the 
ellipticity eszE = 1— A?_ /A+, calculated using the signal be- 
tween the second and third contours, is eszE = 0.20 (0.24) 
for the SCDM (ACDM) cluster. The ellipticities calculated 
using the gas density alone are 0.21 and 0.36, respectively. 
Thus, analytic models to calculate eszE assuming isother- 
mal gas will err by a large margin in the presence of a merger, 
even if eszE is calculated outside the core region. This is un- 
like what we have found for X-ray ellipticities. 

Furthermore, even in the absence of a merger, the 
changing triaxiality of the dark matter halo makes model 
predictions for ellipticity in the SZ maps miss the values 

the solid contours, in order to compare shapes at a given radial 
distance. 



Figure C3. Temperature decrement map for the two clusters of 
Figures IC1I and IC2I in the yz-plane of the corresponding simu- 
lation boxes. The axes are in pixels, and each box is 1.25 h~ 1 
Mpc across. The solid lines show contours of constant J pgas T, 
spaced by factors of 3. The innermost solid contour corresponds 
to a level of 60% (50%) of peak value for the SCDM (ACDM) 
cluster. The top panels compare the shape of the contours of con- 
stant J Pgas (dashed lines) to the decrement-level contours. The 
bottom panels compare the prediction of the analytic model for 
/ Pgas (dotted lines) with the decrement-level contours. See the 
text for further discussion. 



eszE by a larger margin than in the case of X-ray elliptic- 
ity. The bottom panels of Figure IC3I show the predictions 
of the analytic model for J p gas (dotted lines) compared to 
the "isodecrement contours" of the top panels (solid lines). 
It can be seen there that, even in the absence of a merger 
(left), the ellipticity of the analytic model contours is too 
large (even if compared to the simulation-data contours for 
/ pgas only (dashed-lines of top panels). For example, the 
ellipticity between the second and third contours of the ana- 
lytic model predictions is 0.33 (0.46) for the SCDM (ACDM) 
cluster. We find similar results for the other LOS. Thus, re- 
liable predictions (i.e. within 10% of gas-data values) for 
ellipticity in SZE maps need to incorporate the changing 
triaxiality of the dark matter haloes. However, the model is 
still useful to predict statistics of cluster samples such as the 
mean and dispersion. See Section [5] 
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